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Abstract 

The notions of bienergy of a smooth mapping and of biharmonic map 
between Riemannian manifolds are extended to the case when the do- 
main is Finslerian. We determine the first and the second variation of 
the bienergy functional, the equations of Finsler-to-Riemann biharmonic 
maps and some specific examples. We prove that two notable results in 
Riemannian geometry concerning the inexistence of nonharmonic bihar- 
monic maps still hold true this case. 
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1 Introduction 

Biharmonic mappings, as a generalization of harmonic ones, are among the most 
important mappings in physics; initially appearing from problems of elasticity 
theory and fluid mechanics, [24j, in the latter decades, they proved to be useful 
also in computer graphics, geometry processing, |12| and radar imaging, [T]. 
Mathematical arguments, [17:, for the use of biharmonic maps include the fact 
that harmonic maps do not always exist - and biharmonic maps can "succeed 
where harmonic maps have failed"- together with stability issues. On the other 
side, Finslerian models seem to gain more and more ground wherever anisotropy 
of some kind is involved (and not only), from domains such as: kinematics, 
elasticity theory, [6j, seismic ray theory, [3], [29], [30], gravity theories, [TT] . 
[19] . [28] . geometrical optics, [2], thermodynamics, statistical mechanics, [2], 
[H], and up to biology 0, [1]. 

While in Riemannian geometry, biharmonic mappings have been quite in- 
tensively studied (see, for instance, [S], [5], [T7], [32], [33]), to our knowledge, 
in Finsler geometry, only harmonic maps have been considered, [14] . |15| . |16j . 
[T5] , [20] , [15] . Still, the rich potential of Finslerian geometric models makes us 
think that such a study is at least necessary. 

As a first step in this direction, we study in this paper biharmonic mappings 
having as domain real Finsler spaces {M, g) and as codomain, Riemannian ones 



(M, g) or, briefly, Finsler-to-Riemann mappings. Our study continues the one 
made by Xiaohuan Mo and collaborators ([H], [15], [16]) for Finsler-to-Riemann 
harmonic mappings. 

First of all, we extend the concept of bienergy functional for Finsler-to- 
Riemann mappings and determine its Euler-Lagrange equations, i.e., the equa- 
tions of Finsler-to-Riemann biharmonic maps. This process points out a gener- 
alization of the rough Laplacian from Riemannian geometry. 

Any Finsler-to Riemann harmonic map is biharmonic. Just as in Riemannian 
geometry, there exist several cases in which the converse is also true; two notable 
results in Riemannian geometry, due to Guoying Jiang, [5], and C. Oniciuc, |22] . 
respectively, can be generalized without difficulty to our situation: 

1) Any biharmonic mapping whose domain is compact and boundary less and 
whose codomain has nonpositive sectional curvature, is harmonic. 

2) Any biharmonic mapping whose codomain has strictly negative sectional 
curvature, obeying the conditions: a) the norm of its tension is constant and b) 
its rank is greater or equal to 2 at least at one point of its domain, is harmonic. 

Further, we study the biharmonicity of the identity map id : (Af , g) — 
(M, g) in two cases of Finsler-to-Riemann transformations of metrics g t-^ g, 
thus pointing out examples of nonharmonic biharmonic maps. The second case, 
that of linearized perturbations, is inspired from general relativity; even though 
we only consider here positive definite metrics, in our opinion, it is illustrative. 

In the last section, we determine the second variation of the bienergy func- 
tional. Except for the facts that each of the expressions of the tension and of 
the rough Laplacian gains an extra term and of the use of nonlinear connections 
on TM, the first and second variation of the bienergy remain formally similar 
to their Riemannian counterparts. 

In the study of a Finsler space {M, g), there are two major - and equivalent 
- approaches: the one based on the tangent bundle {TM,tt,M), via horizontal 
hfts, and the one based on the puUback bundle tt*TM. As noticed in [2^, the 
study of harmonic maps between real Finsler manifolds is usually carried out on 
Tr*TM (as in [T3], [H]) while in the case of complex Finsler manifolds, it relies 
on the geometry of TM. In order to obtain a more unified method, we preferred 
to work, also in the real case, on TM; the geometric structures we used are the 
TM-correspondcnts of those in [T4]. [T5]. [T6]. 

2 Biharmonic maps in Riemannian geometry 

In this section, we present in brief some results in [T7], [5], 0- 

Let (M, g) and (M, g) be two C°°-smooth, connected Riemannian manifolds 
without boundary, of dimensions n and n; unless elsewhere specified, we will 
assume, as in 0, that M is compact and orientable. On the two manifolds, we 
denote the local coordinates by {x'')i^Yn' (^°')a=Tn^ Levi-Civita connections 
by V, V (with coefficients F'^,^., f"^^) and by r(E), r{E), the modules of C°°- 
smooth sections of any vector bundles E, E over M and M. Commas and 
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will mean partial differentiation with respect to and 5" and 9,;, 9q,, the 
natural bases of the modules T{TM) and T{TM) respectively. 

A C°°-smooth mapping (f) : M ^ M \s called harmonic, if it is a critical 
point of the energy functional 



E : C°°{M, M) M, E{(j)) = 




M 



where d(f) is regarded as a section of the bundle T*M ® cj) ^TM, \\d<j)\\^ = 
traceg{(j)*g) = g"^^ cjapcjf i4>^ j is the squared Hilbert-Schmidt norm of c?(^ and 
dVg is the Riemannian volume element on M. 

Harmonic maps are solutions of the equation t((/)) = 0, where, |17] . 

t{4>) = g'^ {V^dcj^idj) - d0(Va.9,)} g'\V%d4>)dj, (2) 

is a section of the bundle (j)~^TM, called the tension of and is the con- 
nection induced by V in the pullback bundle (f)~^TM, [5'. In local writing: 

r"(0) = g^nr.,, + O,/,'^",, - r%-r, J. 

The above notion of harmonicity generalizes the usual one for mappings be- 
tween Euclidean spaces; notable examples include geodesic curves and minimal 
Riemannian immersions. 

Biharmonic maps (j) £ C°°{M,M) are defined as critical points of the hi- 
energy functional: 

E2{<ly)^\j {T{4>),T{4>))dVg- (3) 

M 

here ( , ) denotes the scalar product on the fibers of TM, determined by g. 
The Euler-Lagrange equation attached to the bienergy is, [TTj : 

r2(0)=O, (4) 

where the bitension T2{4>) of (/> is the section of the bundle (j)~^TM given by: 

T2{<f>) = -A'^T{(t,) - traceg{R^ {d(l),T{<f>))d<f>) (5) 

and the operator A"* ~ ^9^^ ^ti^tj ~ ^Va dj) ^^"^ rough Laplacian, acting 
on sections of (jy^^TM. 

Remarks: 1) Equation ^ is the Riemannian generalization of the bihar- 
monic equation in Euclidean spaces. J24|. 

2) Any harmonic map (p : AI M in biharmonic. 
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3 Finsler structures 



In the following, except for the metric structure on M (and related quantities) 
we preserve the notations and conventions in Section 2. We denote by TM and 
TM the tangent bundles of the manifolds M and M and their local coordi- 
nates, by [x^y] :— (x',y*), (i, y) dots .i and will mean partial 
differentiation with respect to y' and y". 

A. Metric structure: A Finsler structure, [7], [16], on the manifold M is 
a function F : TM -> R with the properties: 

1) F{x,y) is C°°-smooth for y ^ and continuous at y = 0. 

2) F{x,Xy)^XF{x,y),yX>0; 

3) The Finslerian metric tensor: 

g^j{x,y) ■.^^{F^{x,y)).ij (6) 

is positive definite. 

The arc length of a curve c on the Finsler space {M,g) is given, [3, by: 

l{c)= [F{x,dx). (7) 



Condition 2) above insures the independence of l{c) of the chosen parametriza- 
tion of c. 

B. Nonlinear connection and adapted frame on TM : Ehresmann (or 

d 

nonlinear, |7j) connections TTM = HTM ®VTM, with VTM = Span{ — ) 

help simplify computations in Finsler geometry and obtain geometric objects 
with simple transformation rules. A typical choice is the Cartan nonlinear 
connection, built as follows, [?]■ 

Geodesies of the Finsler space (M, g) are defined as critical points c of the 
arc length ([7]); in the natural parametrization, their equations are: 

^ + 2G\x,y)=0, y = x; (8) 

with 2G''-{x,y) = ig''' {{F'^)-h,ky^ - {F'^),h) ] this defines the local coefficients 
G'j- = G^j{x, y) of the Cartan nonlinear connection as: 

G'j-.^G'.j. (9) 
The Cartan nonlinear connection gives rise to the adapted basis: 

(10) 
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on T{TTM) and to its dual {dx\ Sy' = dy' + G'jdxi). 

With respect to coordinate transformations on TM, induced by coordinate 
transformations = (x) on M, the elements of the adapted basis/cobasis 
transform by the same rules as vector/covector fields on M, [5]. 

Any vector field X on TM can be decomposed as: X = hX + vX, hX = 
X'Si, vX = X^di] its horizontal component hX and its vertical component vX 
are vector fields on TM. This leads to a simple rule of transformation for X*, X*. 
A similar situation holds for 1-forms w = huj + vuj, huj = LOidx^, vuj = uJiSy^, [5]) 
and more generally, for tensors of any rank on TM. 

Using Cartan nonlinear connection, tangent vector fields to lifts c' ;= (c, c) 
to TM of unit speed geodesies of M are always horizontal, [7], [S]. 

The adapted basis {6i, di} is generally nonholonomic, the Lie brackets of its 
elements are: 

[Sj,Sk] ^ R'jkix,y)d„ [Sj,dk] = G'jkix,y)d.„ [dj,dk] = 0; 

where: 

K'^ki^,y) - - d,G\, G^,(x,2/) := G\^.^ix,y). (11) 

C. As covariant differentiation rule on TM, we will use the one given by 
the Chern-Rund affine connection D on TM, [6] , locally described by: 

Dsjj = r,kS^, Ds^d, = r^,d,, DgJ, = Dg^d, = 0, (12) 

where = -g^'^{5kghj+Sjghk—Sh9jk) are the " adapted" Christoffel symbolic 
■'2 

of g. The Chern-Rund connection preserves by parallelism the horizontal and 
vertical distributions on TTM, i.e., 

DxihY) = h DxY, DxivY) = v DxY; 

it is generally, only ^.-metrical: 

Dhxg = 0, VX G r{TTM). 

The Chern-Rund connection D has nontrivial torsion: 

T = R'j^d, ® dx^ ® dx^ + P'jkd, ® Sy^ ® dx^ , (13) 

with R^jk as in (ITT]) and P*^^, = G'^j, — T'^j,; the latter defines a horizontal 
1-form: 

P^Pdx', P.-^P%, (14) 

^In the Finslerian case, the usual Christoflfel symbols 7* j. = —g^^ighj,k + 9hk,j ~ 9jk,h) 
do not generally represent the coefficients of an affine connection on TM. 
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which win be used in the following. We notice that the torsion of D has only 
vertical components: 

hT{X, Y) = 0, VX, Y e V{TTM). (15) 

The curvature i? of D is locally described by: 

R = R^\idx'' ® dx^ ® dx^ ® 5^ + Rj\idx'' ® dx'' ® 5y^ ® ^^+■, (16) 
+Pj\iSy^ «) dx'' (g) dx^ (»S^ + P/ki^y^ dx^ ® Sy^ ® di, 

where i?/,, = s^r^, - <5,r*^, + r'^,r\; - r''^,r\, and p/,, = r^^.,,. 

We consider the above notions also for the Riemannian manifold (M , g) and 
designate them by tildes. In this case: f"^^ = G"^^ = ^"p-^, that is: 

[VxYf ^ b^KY'\ yX,Y £r{TM), (17) 

where the superscript indicates the horizontal of vector fields from M to TM; 
also, G°'p{x,y) = ^°'p^{x)y'^ , P"^^ = 0; the Chern-Rund connection becomes 
"fully" metrical: 

Dxg = 0, VX e T{TTM). (18) 
The only nonzero local component of curvature tensor R is i?^"^^, i.e.: 

R{X, Y)Z = R{hX, hY)Z, yX, Y,Z e T{TTM); (19) 

Rp'jS coincide with the components of the curvature R^ of the Levi-Civita 
connection of g and are thus subject to the same symmetries; Bianchi identities, 
[7], also acquire the same form as those of R^ . Ricci identities of D take the 
local form: 

Ds Ds Z'^ - bs bs Rf?.,„Zl^ + R^ioZ"a; , , 

bg^bs^z'^-bs^bg^z'^^o, yz = z''6c,er(HTM). ^ ' 

Another useful property of D is: 

Ds,r = 0. (21) 

The Riemannian metric g gives rise to a scalar product on the fibers on 
HTM, which will be denoted by ( , ) : 

{X, Y) = gapX°'Y^, yX = X°'~Sa,Y = Y^S^ G T{HTM). (22) 

D. Volume form and integration domain. Consider the Riemannian 
volume element dVg = v det Gdx A Sy — det gdx A Sy on TM, determined by 
the Sasaki lift, 0, G := gijdx^ (g) dx^ + gijSy^ (g) 5y^ of g to TM. For an x, y- 
dependent function a : TM — > R, we will consider its integral on the total space 
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of the unit bah bundle of M (which is compact, since g is positive definite), 
divided by the volume VolW^ of the unit ball in the Euclidean space M", i.e., 

a{x, y)dVg = y^j^, j ( j a.{x, y) det g{x, y)dy)dx, (23) 

BM M 

where — {y & T^M \ gijy^y^ < l} (on M, this construction provides a gen- 
eralization of the Riemannian volume element, called the Holmes-Thompson 
volume element, [25]). 

The divergence divX = SiiX"^ det g) — G"' X', of a horizontal vector 

det g 

field X = X'''5i on TM, [21], can be expressed in terms of Chern-Rund covariant 
derivatives as: 

divx = Ds^x' -p{xy, 

we will also use this relation in the form: 

g'^S^Xj = dtvX + g'^r%Xk + P^X\ (24) 

4 Some remarks on Finsler-to-Riemann maps 

Let (j) : M ^ M, (x*) (0"(a;*)) be C°°-smooth. Between the tangent bundles 
TAI and TM, it acts the differential $ := dcj) (regarded as a mapping between 
manifolds); throughout this section, we will use alternatively the two notations 
$ and d(j). The connection D determines a connection D'^'^ in the puUback bundle 
d(l)''\TTM): 

D^^($*y) := D^,xY, yX e r{TTM), Y e T{TTM)- (25) 

with ($*F)(:r,j^) = hcuCC, 

Df{hZ):=:hDfz, Df{vZ)^vDfz, \/Z (^T{d4)-\TTM)). (26) 
The mapping $ is locally described by: 

$: i" = 0"(x), r = '^",(^)y^'. 

With 0"'(a;,y) := ^"^(a;)?/^ we have, for X e T{HTM) : X(0")a„ + 

X{(t)" )da- Taking X = X^6i and expressing in the adapted bases, 

$»X = + (X^5,0"' + iV"^X(/))9„. (27) 

The horizontal component h^^X will have a peculiar importance. 
Lemma 1 For any horizontal vector field X = X^di on TAI : 

m,X = =: d(jJ\X), (28) 

where dcf)^ := 0" dx' 5^ is the horizontal lift of the vector-valued 1-form d(j). 
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Consider a l-parameter variation f : x M, f = f{e, x), /(O, x) — of 
(j) and: 

F -.^df : T{I, X M) TM. 

On T(/£ X M), the local coordinates are {e,x,e' ,y). Taking on the interval 
le C M, the Euclidean metric and the product Finsler metric on x M, we 
will obtain a trivial prolongation of the Cartan nonlinear connection to this 
new manifold, which produces the adapted basis {de,Si,de' ,di} and a trivial 
prolongation of the Chern connection D (which we will denote again by D). i.e., 
Dg^Si = Dg^di = 0, Ds-de = Dg de ~ etc. We also notice that [de^Si] — 

and [dg.di] = 0. 

The connection D'^f will be prolonged to df-^{TTM), by: 

Dgl{F*X) := Dp^g^X, XeT{TTM). (29) 

Lemma 2 For any X,Y e T{TTM), Z e T{df-^{TTM)), p = {x,y) £ TM: 

Dl{df\X)) = Df{df\d,)). (30) 

Proof. (dSl) says that = hf{F^de,F^X) = h{Df^{F^X) ~ D%{F^de) - 
[F^X, F^de]}. The resuh follows then from ^ and [X^de] = 0. ■ 



5 Bienergy and its first variation 

The energy of a Finsler-to-Riemann mapping : Af — M, [14], is i?((?I>) = 
77 / 9^'' 9ai34'%i'P^ jdVg, (note that = (j>{x), hence = Si(j)"), or, in our 
language: 

E{cf>) = i y (d(/.''(5,),d'^''('50) dVg, 

BM 

with ( , ) as in ([2^ . The tension of : M — M, [2], can be regarded as a 
section of (dc/))-'^ HTM: 

t{4>) = g^HD'',fdcf/\S,) - d/'(I?5.<5,) - P^dcf>\d,)}. (31) 

The mapping is harmonic iff its tension vanishes identically. 

It appears as natural to define the bienergy of a mapping </> : M — ^ M as: 

^2(0) = ^ I (T(<^),r(0))dV,. (32) 

Accordingly, by a Finsler-to-Riemann biharmonic map we will mean a critical 
point of the bienergy (f32l) . 
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In order to determine the critical points of E2, we take variations / = /(e, x) 
of (f) as above and denote by 

V := = V:=V\,^o. (33) 

the horizontal lift of the associated deviation vector field df{d^). 

Since the Chern-Rund connection on the codomain TM is metrical, 

^^^'^^ItsJ (-if)^-U))dV,^ J {Dfjif),rif))dV,. 

BM BM 

Let us evaluate the term D'^^r^f) : 

Dtr{f) = Dl{g^={Dl {df\5,)) - df\Ds^5,) - Pdf\5,))} = 
= g^^ [dIdI {dfHs,)) - Dlidf^Ds^S,)) - D^mdfHs,))} 

(g*-' can be taken in front of the Sg-derivative, since in g*-' = g^^{x,y) the 
coordinates x,y do not depend on s). Commuting derivatives by means of the 
curvature tensor of D, taking (fT9)) and [5i, d^] = into account, 

DfDf^ {df\5,)) = i?(V, df\5m''{5,) + DlDl{df\5,)). 

By dSO]), the term Df^Dfjdf~^{dj)) becomes ofofY. Using ^ also in the 
expression D'Q^{df^{Ds-dj)) — D'Q^{Pidf'^{6j)) and summing up, we get: 

= g^HRiV,dfHs^MHSJ)+DfDfV- (34) 
-DtjV-P^DfV}. 

We notice the operators 
g^^-D'^^Df^+D'l ,^+P,D'^^)^:A''f.J = -A''f-tracegR{df\ ■ )df\ (35) 
acting on sections of the bundle {df)^^{HTM). With this, we have: 

dItU) = -A'^-^V - g^=R{df\5.l\)df\5,) = ^^(V). (36) 
Evaluating at e = and substituting into the expression of the variation, 

^(/)le=o= / {J{VlT{ct>))dV,. (37) 

BM 

It remains to transform the above expression so as to have V in the right hand 
side of the scalar product. This will be easy using the following lemma. 
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Lemma 3 The operators A*^"^ and J are self-adjoint: 



{A''*X,Y)dVg^ I {X,A'''^Y)dVg, I {jX,Y)dVg^ I {X,JY)dVg 

(38) 



BM BM BM BM 



for any X,Y e T{d(t>^^{HTM)). 

Proof. We start from the left hand side of the first relation (1551) ; integrating 
by parts the term J (— g^^ of, of, X , Y) dVg and applying (p4)) . we get: 



BM 



(A'^'l'X, Y) dVg = - J g'^ (D'^fx, D'^fY^ dVg. (39) 

BM BM 

Integrating once again by parts, we obtain The self-adjointness of J 

follows then from the symmetries of i?. ■ 

The operator A'*'^ is a generalization of the rough Laplacian from Rieman- 
nian geometry, built in the same spirit as the horizontal Laplacian acting on 
differential forms in [31], [32] . 

Using Lemma [3] in (|37l) . we get: 

Proposition A a) The first variation of the bienergy of a mapping (j) : A'l ~> M 
from the Finsler space (M,g) to the Riemann space {M , g) is: 

^^^|.=o= J {-A^*T{^)-tracegR{dcjJ',Tmd<fJ',v)dVg- (40) 

BM 

b) The mapping (j) is biharmonic iff: 

r2(0) := -A'*'^r((/.) - tracegR{d<t)^ ,T{4)))d(t/' = 0. (41) 

Remarks. 1) In the above, we considered, as in ^8 , that M is compact 
and without boundary. Elsewhere, all the discussion can be made on a compact 
subset T) of M] in this case, we assume that, on the boundary of I?, the vector 
field V and the covariant derivatives Df^V vanish. 

2) Any harmonic map from a Finsler space to a Riemann one is biharmonic, 
namely, a minimum point for the bienergy functional. A biharmonic map which 
is not harmonic will be called proper biharmonic. 

Particular cases: 

1) If M = M" with the Euclidean metric, then the biharmonic equation (|4ip 
becomes: 

A'^'^r(0) = 0. 

2) If M = S" is the unit Euclidean sphere, then, using the expression of the 
Riemann tensor of a space form, we get that (j) : M ^ §" is biharmonic iff: 

A'^'t'TicI)) + 2e{(l))T{(f>) - traccg (d(f>'' ,t {(/>)) dc/)'' = 0, 
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where e((/)) — -traccg (^d(j)^ , dcp'^J is the energy density of (p. The resuh is similar 

to the one in the Riemannian case, [5]. 

3) If M is a weakly Landsberg manifold, i.e., [Ij, if P = 0, then the ex- 
pressions of the tension and of the rough Laplacian become formally similar to 
the ones in the Riemannian case - just, depending on the fiber coordinates y*: 

6 Existence of proper biharmonic maps 

The following two results represent generalizations to Finsler-to-Riemann maps 
of two theorems in [8 and 22 respectively. 

Theorem 5 // (M, g) is a compact Finslerian manifold without boundary and 
{M,g) is Riemannian with nonpositive sectional curvature, then any biharmonic 
map (j) : M ^ M is harmonic. 

Proof. The proof follows similar steps to the one in the Riemannian case, [5]. 
We apply the horizontal Laplace-Beltrami operator, |31| . A/ := —div{gradh /), 
where gradhf '■= {g^-'5jf)5i, to the scalar function / := ||r((/))||^, defined on 
TM : 

-iA ||r(0)f = \{Ds^{g^^5, ||r(0)||^) - g^^P,6, ||r(0)f } = 
= \9H6^5, hm? - Wrmf - P,S, ||r(0)f }. 

Taking into account that ^\jSk = Dg.Sj and expressing the action of the 
adapted basis vector fields 5i,6j on ||t((/))||^ = {t{(P),t{4>)) in terms of D'^'^- 
covariant derivatives, we obtain: 

- Ia \\rm' = - (A'^M^), r(0)) + g^= (^Di'^T{^),Di'^r{cp)) . 

By means of the biharmonic equation (1411) . this becomes: 

- \A\\Tm^ - ltrace,R{d(p'\T{4>))d^^.r{(p))+g^^ (^D'',K{4>),d1K{4>)) . 

/42) 

According to the hypothesis that the sectional curvature of (M, g) is nonpositive, 
the curvature term above is nonnegative; since g^^ (^Dg'^T{(j)), Dg'^Tlc/))'^ (as a 

1 2 

squared norm) is nonnegative, too, we get: ~ 2^ ll'''(<^)ll — 0- 

On the other side, we have, EB, / A ||T((/.)|r dVg = 0, hence, A ||t(0)|P = 

BM 

0; thus, by (|^ . .g'-' (^DgfT{(j)),Dg'^T{(j))'^ = 0; as a consequence, 

Dtfri^) - 0. (43) 
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Take the horizontal vector field X {g^^ ((^ j, t((^)))5j on TM; by (gS]), we get: 

0= y" divXdVg^ j {T{(t>),T{(t>))dVg 
BM BM >o 

and therefore, {t{(I)),t{4i)) = t(0) ~ 0, i.e., cj) is harmonic. ■ 

Dropping any condition upon the compactness or on the boundary of Af , we 
have: 

Theorem 6 Let {M,g) be an arbitrary Finsler space (not necessarily compact), 
{M,g), a Riemannian manifold with strictly negative sectional curvature and 
(f) : M ^ M, a biharmonic map. If (j) has the properties: 1) ||t(0)|| — const, 
and 2) there exists a point xq G M at which the rank of (j) is at least 2, then (j) 
is harmonic. 

Proof. The proof is similar to the one in the Riemannian case, [12] • From the 
hypothesis ||t((/))|| = const., in (|42|). the left hand side is 0; but both terms in 

the right hand side are nonnegative, hence: (trace g R{d(j/^ , T{(j)))d(j/^ , = 0. 

Since Rierug < 0, we must have, for all i = l,n : d(f)^(5i) \\ t{(I)). Taking into 
account that at a;o, rank{(j}) > 2, the only possibility is T{(j)){xo) — 0. Using 
||t((/))|| = const, it follows that t{4>) = 0, i.e., 4> is harmonic. ■ 



7 Biharmonicity of the identity map 

Throughout this section, we assume that M = M (not necessarily compact), 
dimM = n, and denote the coordinates on TAI by {x'^,y^). Considering on M 
two metrics: a Riemannian one g and a Finslerian one g, we will explore the 
biharmonicity of the Finsler-to-Riemann mapping: 

id:{M,g)^{M,g). (44) 

In this situation, there appear two adapted bases {5i, di) and {5i, di) on TM, 
together with the covariant differentiations given by D, b and According 
to [16j , the tension of the identity map has the local components 

r\^d)^g^\^^^^G^^^) (45) 

(note: our is half the one in [16]). 
Let us denote b -.^ — F"^ , i.e.: 

9i] y) = gij {x) + b,j {x, y), (46) 

where the function b = b[x, y) is homogeneous of degree 2 in y and 6^ ~ b.ij. 

In the geodesic equations ([8]), we express the derivatives in terms of Dj - 
covariant derivatives, denoted in the following by double bars ||i ; we obtain: 

2G' = 2& + 2B\ (47) 
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where: 



1 

2.y v^y/i||jf " ~ \\h 

1 



2B^ -.^ -g^\2v,,\\,y' ~ F\,:) (48) 



and yh '■— l^P^h — dhjV'' ■ The tension of id is: 

T^{id) = (49) 

A direct computation shows that, in (1481) . the covariant derivative 2i/^||j can 
be rewritten as: 

'^yiA\j=F\rh- (50) 
Remarks: 1) If b is parahel with respect to Z?, then + ~ 0? 

taking into account ((50|) . we get _B* = r* = 0; in this case, the identity map 

is harmonic, i.e., also biharmonic. 

2) Assuming that g is a Berwald-type metric, i.e., G'^j, = G'*jj,(x), then 

there exists, [13], [27], a Riemannian metric such that G'^j. = G*^^.; thus, the 
identity map is, again, harmonic, hence, biharmonic. 

We win find in the fohowing two examples of Finslerian perturbations h for 
which the identity of M is proper biharmonic. 

With T* := T^{id), the relation between the and Z)-covariant deriva- 

tives of is: 

dX^t^ = 5,T^ + f V,r^ = r^ii^ - B\t\^. (51) 
1. Suppose that the Finslerian function satisfies: 

F\h^ {a,y)gyh, (52) 

where (a, y) g := gija^y^ and a' = a*(a;) are components of a vector field A — a^di 
on M (relations ((52]) are equivalent to a first order ODE system in g^ ). 

A brief calculation using (1501) leads to: 2yj^^hy^ ^ -P'^ij = o.jF'^i that is, 
2S* = a*(j:)F2. From gH]), we obtain: 

T' = -]^na\ (53) 



Since r* = T'(a;), relation (|5T]) becomes simply: D r* = r*||^ and the 



biharmonic equation is written as 

5''(^^imi/c-rVii/-^//fc^') = o. (54) 

Here, taking into account that Rjuk — Rikij and Ricci identities (|20]) for £), the 
curvature term Rj i}.t^ can be expressed by commuting i'-covariant derivatives 
of t\ It turns out that a sufficient condition for the biharmonicity of id is: 

r^ll^. = 0. (55) 

{Note: this statement is always true in the Riemannian case, but generally, 
not in the Finsler-to-Riemann one, where, as a rule, r* = T^{x,y)). 
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Using ([55]) . we deduce that ([55]) is identically satisfied if the vector field 
A'^ = a^Si is parallel with respect to D. But, according to ([T7|) . this is nothing 
but: Vd^A = 0. In other words: 

Proposition 7 //, in I152\) . the nonzero vector field A = a^{x)di is parallel with 
respect to g, then the identity map id : (M, g) — >■ (M, g) is proper biharmonic. 

2. Linearized Finslerian perturbations of the Euclidean metric. 

Assume that (M,gij) = (M",(5y) and the perturbation bij =: eij{x,y) is small 
(linearly approximable) , that is, we may neglect all terms of degree greater than 
one in Sij and its derivatives, |10j . In this case, the inverse metric is given by: 
gik ^ jjfc _ ^ik g^j^j relation (gS]) becomes: 

We notice that the tension r will be of the same order of smallness as e; 
it means that products of r with e and its derivatives can be neglected. For 
instance, we have: B^t"^ ~ 0, which, substituted into (ISD) . leads to: 

The biharmonic equation takes the simple form: j'™r' / „ = 0. Again, a suffi- 
cient condition for biharmonicity is 

(or: = r'(j/)), that is, S''' {shj,k,i + £hk,j,i - £jk.h,i)y^y^ = 0. We obtain: 

Proposition 8 Let the Finsler metric gij{x,y) — 6ij + eij{x,y) be a linearized 
perturbation of the Euclidean metric on R". If the components Sij{x,y) are non- 
constant and affine in x, then id : {M",Sij) — > (IR",^) is proper biharmonic. 

8 Second variation of the bienergy 

Take a biharmonic map (j) : (Af, g) — >■ (M, g) and a smooth 2-parameter variation 
/ = /(ei, £2, x), /(0, 0,x) = (l) of (/), with 

V.^dfHde,), V2=d/(aeJ, := Vi|,,=,,=o, V^2 := Vi|,,^,,=o 

(if M has a boundary, then Vi , V2 and their ^^-covariant derivatives are assumed 
to vanish on dM). 

The deduction of the second variation of E2 follows the same steps as in 
the Riemannian case, with two differences: in the expressions of t(/) and of 
A''-'', there appear extra terms and we have to take into account that $*(9(;J 
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is, generally, not horizontal. Fortunately, as we will see below, these will finally 
not complicate the expression of the variation. 

We denote, for simplicity, r := t(/). According to ([iD|) . (|1T|) : 



BM 



J (r2(/),Vi)dV,; (56) 



differentiating with respect to £2 



deide2 

BM 

At Si = 62 = 0, since (f> is biharmonic, the second term in the right hand side 
term will vanish. It is thus enough to evaluate the first one; we have: 

Dtj2{f) = -Dfj^'^r) - D%^{trace,R{df\T)dh. (57) 

The covariant derivative of the Laplacian — A'^-'^t is: 

Ti := -Dl^ (A'^fr) = g^Wf ^ {d%D%t- d\,t P^Dfr)) . (58) 



Commuting covariant derivatives by means of the curvature R (twice for the 

V r)df , 

Si Sj ' 



term D'^g D'^g t), taking into account that [de,Si] — and ([T5)) . we find: 



The terms in D'g r can be grouped into —A'^f {D'q t) : 

Ti = ~A''f{Dlj)+g^={R{\2,df\5,))D%j+ 
(^R{V2,df~HSj))T) - R{V2,df^{Ds,6,))T - P,R{^2,df^{6,))T}. 

Splitting D'^g, (^R(V2,df'^{Sj))Tj as a sum of derivatives, we recognize in the 
resulting expression i?(V2, t)t : 



Ti = -A'^fiDilj) + R{Y2,t)t + g^^{iD%R)iY2,dfHs,))T- 
+2R{V2,df~\S,))D%T + R{D%Y2,df'\S,))T}. 



(59) 



The curvature term T2 := ~Df^^{tracegR{df~^,T)df'') in dST]) is: 

T2 = -^g'^iDf^^R){dfHs^),T)dfH6,)+RiDi{^dfHs^),T)d^^^ 
+R{df\S.),Dl T)df\5,) + R{df\5.),T)D% (dfHSj)). 



(60) 
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Taking into account that R = R(x) only, we obtain D'q R = D^^^R. Trans- 
forming D'g df^{di), D'q df'^{5j) by ([50)1 and then using first Bianchi identity 
in the second term: 

T2 = -^g'^{{b^Mdf''{h),T)df^{5,) + R{DfY2,T)df\5,) + 

+R{df\5,),Dl^r)df\5,)+R{df>^{5,lT)bs^\2} = 

= ~g^i{{b^,R){df'-{5,),T)df>^{5,)+2R{df^{5,),T)bs,^2- 
~Ridf^{S,),DfV,)T) + R{df^{5,),D%j)df^{5,)}. 

Second, and then first Bianchi identities for the (_Dv2^)"term tell us that: 

-g'^ (Dv, R) (df' {5,), T)df^^ {6, ) = g'HibrR) (V2 , df^ {SiM'^ {5,)- 
~iDs,R){df''iS,),T)Y2 - {bsM^2.df^{5,))T). 

Substituting into T2 and adding: Ti + T2 = D'q T2(/), we get: 



Mf)^ J{D%^ r) + i?(V2 , r)r + g^^ { (brR) ( V2 , d/" (,5,))^/'' ) 
- (7^% R) {df^ {5, ) , r ) V2 + 2i?(V2 , df^ {S,))D%^ r - 2R{df^ {S, ) , T)i^% 

with J' as in (1551) . Using ([55]) and evaluating at £1 = £2 = 0, we get: 



Proposition 9 The second variation of the bienergy of a Finsler-to-Riemann 
biharmonic map 4> : M ^ M is: 



BM 



J {Vi, J^V2 + R{V2,t)t+ (62) 



+g'^{{brR){V2,df'^i6.0)df'\S,) - {D%^R){df\5,),T)Y2 + 
^2R{V2,df''{d,))D%T - 2R{df~''i6,),T)D%V2} ) dVg. 



In particular cases (for instance, M — R" or 5"), (p^ becomes considerably 
simpler. 

d^E2(f) 

The Hessian H : (1/1,^2) ^ n{VuV2) = ^-4^|ei=e.=o of the bienergy is 

C'£iC'£2 

a symmetric bilinear form. A solution of the biharmonic equation is stable if 
the quadratic form 'H(y, V) is nonnegative for any V. As an example, harmonic 
maps are stable biharmonic maps. 
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